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NOMENCLATURE 


a 

— 

crack  depth 

ao 

- 

Initial  crack  depth 

N 

- 

number  of  fatigue  cycles 

C  ,ra 

- 

material  constants  in  crack  growth  equation 

K 

- 

stress  intensity  factor 

kc 

- 

critical  value  of  K 

^max 

- 

maximum  value  of  K  during  fatigue  cycle 

KP 

- 

K  due  to  pressure 

KR 

- 

K  due  to  residual  stress 

AK 

- 

range  of  K  during  fatigue  cycle 

at 

- 

tangential  stress 

atR 

- 

residual  ot 

SY 

- 

yield  stress 

R 

- 

radius 

RO 

- 

outer  radius 

RI 

- 

inner  radius 

RP 

- 

radius  of  elastic-plastic  Interface 

*a 

- 

radius  at  crack  tip 

SA 

- 

°tR  at  R  =  Ra 

SB 

- 

at£  at  R  =  RI  corrected  for  Bauschinger  Effect 

RZ 

- 

radius  at  which  atR  =  0 

SF 

- 

crack  shape  factor 

F 

- 

portion  of  negative  KR  included  in  AK 

Y 

- 

diameter  ratio,  RO/RI 

ill 


INTRODUCTION 


In  many  cases,  the  fatigue  life  of  a  structural  component  can  be 
predicted  with  reasonable  accuracy  by  determining  the  maximum  range  of 
operating  cyclic  stress  intensity;  and  using  this  with  stress  vs.  cycles  (S/N) 
curve,  obtained  from  an  axial  fatigue  test  on  small  specimens,  determine 
fatigue  life.  Fatigue  lives  predicted  by  this  method  for  thick-walled 
cylinders  have  generally  not  been  in  reasonable  agreement  with  the 
experimental  results  (refs  1,2).  This  is  particularly  true  if  the  cylinders 
contain  residual  stresses  produced  by  autof rettage . 

This  disagreement  is  due  to  the  fact  that  a  large  portion  of  the  fatigue 
life  of  a  thick-walled  cylinder  may  be  consumed  by  propagation  of  the  fatigue 
crack.  Another  reason  is  the  fact  that  the  internal  pressure  acts  directly  on 
the  surfaces  of  microscopic  voids  and  defects  on  the  bore  surface  which  may 
decrease  the  number  of  cycles  for  crack  initiation  and  early  growth. 

Based  on  the  above,  it  is  proposed  that  a  more  accurate  prediction  of 
fatigue  lives  of  pressurized,  thick-walled  cylinders  can  be  obtained  from  a 
fatigue  crack  propagation  model  rather  than  a  model  based  on  the  S/N  curve 
approach.  This  has  the  added  advantage  of  being  able  to  readily  account  for 
the  presence  of  residual  stresses  due  to  autofrettage  or  shrink  fitting  or 
other  processes. 


^Eisenstadt,  R. ,  Kendall,  D.  P. ,  and  Davidson,  T.  E.,  "A  Comparison  of  Results 
of  Axial-Tension,  Rotating-Beam,  and  Pressurized  Cylinder  Fatigue  Tests,” 
Experimental  Mechanics,  Vol.  9,  No.  6,  June  1969,  pp.  250-254. 

^Karl,  E.,  "The  Fatigue  Strength  of  Thick-Walled  Cylinders,"  Proceedings  of 
the  9th  AIRAPT  High  Pressure  Conference,  Elsevier,  1984. 
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THEORY 


It  has  been  shown  that  the  rate  of  propagation  of  a  fatigue  crack  can  be 
determined  from  the  applied  range  of  stress  intensity  factor,  ( AK) 9  associated 
with  the  crack  and  the  cyclic  loading  applied  during  the  fatigue  cycle,  by  the 
following  equation. 

da/dN  =  C( AK)ra  (1) 

where  a  is  the  crack  depth,  N  the  number  of  fatigue  cycles,  and  C  and  m  are 
material  constants. 

The  fatigue  life  of  a  thick-walled  cylinder  can  therefore  be  predicted  by 
integration  of  Eq.  (1)  if  the  following  assumptions  are  made: 

a.  The  constants  C  and  m  are  known  for  the  cylinder  material. 

b.  All  materials  contain  crack-like  flaws  which  can  be  assumed  to  be 
initial  cracks  from  which  fatigue  cracks  will  propagate. 

c.  Failure  of  the  cylinder  will  occur  when,  or  very  soon  after,  the 
maximum  value  of  K  exceeds  some  critical  value,  Kc. 

d.  The  range  of  stress  intensity  factor  during  the  fatigue  cycle  can  be 
calculated  as  a  function  of  pressure,  crack  size,  and  residual  stress 
distribution. 

Equation  (1)  can  then  be  integrated  from  an  assumed  initial  crack  depth 
to  the  crack  depth  at  which  K^x  =  Kc,  to  determine  the  fatigue  life. 

RESIDUAL  STRESS  DISTRIBUTION 

In  calculating  the  effect  of  residual  stress  on  stress  intensity  factor, 
only  the  tangential  component  need  to  be  considered  since  the  longitudinal  and 
radial  stress  components  do  not  contribute  to  the  K  value  for  a  crack,  whose 
plane  is  normal  to  the  tangential  direction. 
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Based  on  the  Tresca  yield  criteria  arid  assuming  ideal  elas tic-perfect ly- 
plastic  material  behavior,  the  tangential  residual  stress  distribution  in  an 


autof rettaged  thick-walled  cylinder  is  given  by  (ref  3): 


°tR  RO2  +  RP2  R  RI2  RO2  -  RP2  RP  RO2 

- = -  +-  £n(— ) - 5 - ?  ( - j - +  In  —  )(1  + 

SY  2(R0)Z  RP  R0Z  -  RIZ  2R0Z  RI  Rz 


(2) 


Equation  (2)  is  based  on  the  assumption  that  the  material  behaves 
elastically  on  the  release  of  the  autof rettage  pressure.  However,  most 
materials,  particularly  the  quenched  and  tempered,  low  alloy  steels,  generally 
used  for  high  pressure  vessels,  exhibit  a  significant  Bauschinger  effect. 

This  is  a  reduction  in  the  compressive  yield  strength  resulting  from  prior 
plastic  deformation  in  tension.  Thus  the  material  near  the  bore  surface 
generally  does  not  behave  elastically  during  unloading,  but  actually  "reverse 
yields”  in  compression  which  signif icantly  reduces  the  residual  tangential 
stresses  near  the  bore  from  those  calculated  from  Eq .  (2).  The  effects  of 
this  phenomenon  on  the  final  residual  stress  distribution  has  been  considered 
by  several  investigators  (refs  4-7).  A  general  analytical  approach  to  the 


■^Timoshenko,  S.,  Strength  of  Materials,  Part  II,  Van  Nostrand,  1956,  pp. 
386-392. 

^Franklin,  G.  J.  and  Morrison,  J .  L.  M.,  "Autof ret tage  of  Cylinders: 
Prediction  of  Pressure/External  Expansion  Curves  and  Calculation  of  Residual 
Stresses,"  Institution  of  Mechanical  Engineers,  London,  Vol.  174,  1960,  p. 
947  . 

^Sidebot tom,  0.  M.,  Chu,  S.  C.,  and  Lamba,  H.  S.,  "Unloading  of  Thick-Walled 
Cylinders  That  Have  Been  Plastically  Deformed,"  Experimental  Mechanics,  Vol. 
16,  No.  12,  December  1976,  pp.  454-460. 

^Milligan,  R.  V.,  Koo,  W.  H. ,  and  Davidson,  T.  E. ,  "The  Bauschinger  Effect  in 
a  High  Strength  Steel,"  ASME  Publication  No.  65-MET-9,  1965. 

^Kendall,  D.  P. ,  "The  Effect  of  Material  Removal  on  the  Strength  of 
Autof rettaged  Cylinders,"  Watervliet  Arsenal  Technical  Report  WVT-7003, 
January  1970  (available  from  NTIS) . 
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calculation  of  resulting  residual  stresses  based  on  realistic  modelling  of  the 
actual  material  behavior  has  been  proposed  by  Chen  (ref  8) • 

From  an  evaluation  of  Chen’s  results,  along  with  the  experimental  results 
of  References  4  through  7 ,  the  following  general  conclusions  regarding  the 
actual  distribution  of  residual  tangential  stress  can  be  reached. 

a.  If  the  absolute  value  of  residual  stress  at  the  bore,  calculated  from 
Eq.  (2),  does  not  exceed  0.4  times  the  original  yield  strength,  results  of  Eq . 
(2)  can  be  used  with  no  correction. 

b.  The  distribution  of  residual  tangential  stress  can  be  assumed  to  be 
linear  between  the  bore  radius  and  the  radius  at  which  Eq .  (2)  gives  a  value 
of  -  0. 

c.  The  actual  residual  stresses  at  the  bore  are  significantly  less  than 
those  that  would  be  predicted  from  Eq .  (2).  However,  due  to  the  strain¬ 
hardening  that  is  associated  with  the  reduced  compressive  yield  strength,  the 
actual  residual  stresses  are  related  to  the  calculated  values.  In  other 
words,  the  more  the  calculated  residual  stress  exceeds  the  reduced  compressive 


^Franklin,  G.  J.  and  Morrison,  J.  L.  M. ,  "Autofret tage  of  Cylinders: 
Prediction  of  Pressure/External  Expansion  Curves  and  Calculation  of  Residual 
Stresses,"  Institution  of  Mechanical  Engineers,  London,  Vol.  174,  1960,  p. 
947  . 

^Sidebottom,  0.  M. ,  Chu,  S.  C.,  and  Lamba,  H.  S.,  "Unloading  of  Thick-Walled 
Cylinders  That  Have  Been  Plastically  Deformed,"  Experimental  Mechanics,  Vol. 
16,  No.  12,  December  1976,  pp.  454-460. 

^Milligan,  R.  V.,  Koo,  W.  H. ,  and  Davidson,  T.  E. ,  "The  Bauschinger  Effect  in 
a  High  Strength  Steel,"  ASME  Publication  No.  65-MET-9,  1965. 

^Kendall,  D.  P.,  "The  Effect  of  Material  Removal  on  the  Strength  of 
Autof rettaged  Cylinders,"  Watervliet  Arsenal  Technical  Report  WVT-7003, 
January  1970  (available  from  NTIS) . 

^Chen,  P.  C.  T.,  "Prediction  of  Residual  Stresses  in  an  Autof rettaged 
Cylinder,"  Proceedings  of  the  9th  AIRAPT  High  Pressure  Conference,  Elsevier, 
1984. 
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yield  strength,  the  more  reverse  yielding  or  plastic  strain  will  occur  on 
unloading.  This  will  result  in  an  increase  in  compressive  yield  strength  due 
to  strain-hardening. 

Based  on  the  above  observations,  the  following  equation  is  proposed  for 
the  residual  tangential  stress  at  the  bore,  corrected  for  the  Bauschinger 
effect,  defined  as  (SB). 


SB  =  for  atRI  >  -0.4  SY 

(3) 

SB  =  -0 . 4 ( SY )  +  0.5  [atRI  +  0.4(SY)]  for  atRI  <  -0.4  SY 
where  o^rx  is  the  value  of  Oj.R  at  R  =  RI  from  Eq.  (2).  The  residual  bore 

t 

stresses  from  Eq.  (3)  are  compared  with  those  from  Chen's  analysis  (ref  8)  in 
Figure  1. 

Defining  RZ  as  the  radius  at  which  the  value  of  atR  calculated  from  Eq. 
(2)  is  zero,  and  assuming  a  linear  distribution  of  o^R  between  RI  and  RZ 
results  in  the  following  equation: 


RZ  -  R 


atR  =  SB( 


- ) 

RZ  -  RI 


(4) 


STRESS  INTENSITY  FACTORS 

A  number  of  investigators  have  published  solutions  for  the  stress 
intensity  factor  associated  with  a  crack  in  a  pressurized  thick-walled 


®Chen,  P.  C.  T.,  "Prediction  of  Residual  Stresses  in  an  Autof rettaged 
Cylinder,"  Proceedings  of  the  9th  AIRAPT  High  Pressure  Conference,  Elsevier, 
1984. 
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cylinder  (refs  9-12).  Much  of  this  work  is  related  to  the  effects  of  residual 
stresses,  crack  shape,  and  multiple  cracks,  all  of  which  significantly  affect 
the  value  of  K.  A  good  summary  of  this  work  can  be  found  in  Parker,  et  al 
(ref  13).  This  reference  also  gives  one  of  the  most  accurate  solutions 
available  for  the  effect  of  autofrettage  residual  stresses  on  straight- 
fronted,  single,  internal  cracks.  Underwood  and  Throop  (ref  14)  previously 
published  an  approximate  solution  assuming  a  linear  varying  residual  stress 
over  the  crack  depth.  For  the  case  given  as  an  example  in  Parker,  et  al  (ref 
13),  the  values  of  K  calculated  using  the  Underwood  and  Throop  solution  agree 
with  the  more  rigorous  analysis  within  a  few  percent. 

Underwood  and  Throop  (ref  14)  give  the  following  equations  for  the  stress 
intensity  factors  due  to  both  the  internal  pressure  (KP)  and  to  the  auto¬ 
frettage  residual  stress  (K.R) 


2.24  PY2  , _ 

KP  =  [-y2 - --](SF)/rra  (5) 

KR  =  [1.12  (SB)  -  0.68  (SB  -  SA)](SF)/ira  (6) 


^Bowie,  0.  L.  and  Freese,  C.  E.,  "Elastic  Analysis  for  a  Radial  Crack  in  a 
Circular  Ring,”  Engineering  Fracture  Mechanics,  Vol.  4,  No.  2,  1972,  pp. 
315-321. 

10Pu,  S.  L.  and  Hussain,  M.  A.,  ’’Stress  Intensity  Factors  For  Radial  Cracks  In 
a  Partially  Autof rettaged  Thick-Walled  Cylinder,”  ASTM  STP  791,  1983,  pp. 
1-194-1-213. 

H-Tan,  C.  L.  and  Fenner,  R.  T.  ,  ’’Stress  Intensity  Factors  For  Serai-Elliptical 
Cracks  in  Pressurized  Cylinders  Using  Boundary  Integral  Equation  Method,” 
Intj^  Jour,  of  Fracture,  Vol.  16,  No.  3,  1980,  pp.  233-245. 

l^Atulri,  N.  and  Kathiresan,  K.  ,  ”3-D  Analysis  of  Surface  Flaws  in  a  Thick- 
Walled  Reactor  Pressure  Vessel  Using  Displacement-Hybrid  Finite  Element 
Method,”  Nuclear  Engineering  and  Design,  Vol.  51,  pp.  163-176. 

l^Parker,  A.  P. ,  Underwood,  J.  H. ,  Throop,  J.  F.,  and  Andrasic,  C.  P.,  "Stress 
Intensity  and  Fatigue  Crack  Growth  in  a  Pressurized,  Autof rettaged  Thick 
Cylinder,”  ASTM  STP  791,  pp.  I-216-I-237. 

^Underwood,  J,  H,  and  Throop,  J.  F.,  ’’Surface  Crack  K-Estiraates  and  Fatigue 
Life  Calculations  in  Cannon  Tubes,”  ASTM  STP  687,  1979,  pp.  195-210. 
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where  P 


internal  pressure 


Y  =  diameter  ratio 

SB  =  residual  tangential  stress  at  bore 

SA  =  residual  tangential  stress  at  the  radius  located  at  the  crack  tip 

a  -  crack  depth 

SF  =  crack  shape  factor  (to  be  discussed  later) 

CRACK  SHAPE 

The  equations  for  the  stress  intensity  factors  (Eqs.  (5)  and  (6)),  are 
for  a  straight-fronted  crack  (constant  depth  over  the  whole  length  of  the 
cylinder).  In  practice  cracks  at  the  bore  of  a  pressurized  cylinder  usually 
grow  in  a  semi-elliptical  shape.  This  shape  will  vary  from  a  semi-circular 
shape  (depth-to-length  ratio  of  about  0.5)  to  a  fairly  long  semi-ellipse  with 
a  depth-to-length  ratio  of  0.1.  Due  to  the  support  provided  by  the  uncracked 
material  beyond  the  ends  of  the  crack,  there  is  a  significant  decrease  in  the 
stress  intensity  factor  with  increasing  values  of  depth-to-length  ratio. 

A  number  of  investigators  have  developed  solutions  for  crack  shape 
factors  for  cracked  plates,  notably  that  of  Newman  and  Raju  (ref  15).  These 
have  been  used  by  Parker,  et  al  and  others,  to  estimate  the  reduction  in 
stress  intensity  factor  due  to  crack  shape.  Tan  and  Fenner  (ref  11)  and 
Atulri  and  Kathiresan  (ref  12)  have  obtained  shape  factors  for  specific  cases 


1^-Tan,  C.  L.  and  Fenner,  R.  T.,  MStress  Intensity  Factors  For  Semi-Elliptical 
Cracks  in  Pressurized  Cylinders  Using  Boundary  Integral  Equation  Method," 
Int.  Jour,  of  Fracture,  Vol.  16,  No.  3,  1980,  pp.  233-245. 

^Atulri,  n.  and  Kathiresan,  K. ,  "3-D  Analysis  of  Surface  Flaws  in  a  Thick- 
Walled  Reactor  Pressure  Vessel  Using  Displacement-Hybrid  Finite  Element 
Method,"  Nuclear  Engineering  and  Design,  Vol.  51,  pp.  163-176. 

■^Newman,  J.  c.  and  Raju,  I.  S.,  "Analysis  of  Surface  Cracks  in  Finite  Plates 
Under  Tension  and  Bending  Loads,"  NASA,  TN1578,  1979. 
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of  cracked  thick-walled  cylinders  and  have  shown  less  reduction  in  K  due  to 
crack  shape  than  is  estimated  from  the  plate  solutions,  particularly  for  deep 
cracks*  These  results  are  discussed  in  detail  by  Parker,  et  al  (ref  13),  and 
show  the  complex  way  in  which  shape  factors  vary  with  crack  depth  and  shape. 

Only  relatively  shallow  cracks  need  be  considered,  since  most  of  the 
fatigue  life  is  in  the  shallow  crack  portion  of  the  crack  growth.  Most 
fatigue  cracks  in  thick-walled  cylinders  have  been  found  to  have  a  depth-to- 
length  ratio  between  0.4  and  0,3,  unless  some  unusual  surface  damage  or 
defects  force  a  lower  ratio.  In  order  to  obtain  a  simple  method  of  calculat¬ 
ing  fatigue  life,  it  is  proposed  that  a  constant  value  of  shape  factor  of  0.7 
be  used  unless  there  is  reason  to  expect  unusually  long  initial  cracks.  In 
this  case,  a  value  of  0.8  or  0.9  can  be  used  for  added  safety.  The  value  of 
0.7  is  probably  conservative  since  actual  values  may  be  lower  than  this  in 
many  cases.  This  will  be  further  discussed  under  Numerical  Results.  This 
shape  factor  is  then  multiplied  times  the  value  of  K  as  shown  in  Eqs.  (5)  and 
(6). 

STRESS  INTENSITY  RANGE 

The  crack  growth  rate,  as  calculated  from  Eq.  (1),  is  a  function  of  the 
range  of  stress  intensity  factor  or  the  difference  between  the  maximum  and 
minimum  value  of  K  during  the  fatigue  cycle  (AK).  If  the  minimum  value  of  K 
is  zero  or  greater  than  zero,  there  is  no  problem  in  defining  AK.  However,  in 

-^Parker,  A.  P. ,  Underwood,  J.  H. ,  Throop,  J.  F.,  and  Andrasic,  C.  P.,  "Stress 
Intensity  and  Fatigue  Crack  Growth  in  a  Pressurized,  Autof ret taged  Thick 
Cylinder,"  ASTM  STP  791,  pp.  I-216-I-237. 
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an  autof rettaged  cylinder,  the  value  of  K  calculated  from  Eq.  (6)  is  negative 
for  shallow  cracks  when  P  equals  zero.  It  is  generally  assumed  that  when  K 
equals  zero  the  surfaces  of  the  fatigue  crack  are  in  contact.  If  a 
compressive  stress  is  applied,  the  crack  surfaces  can  not  close  any  further 
and  thus  the  K  can  not  become  negative.  It  is  common  practice  to  consider 
that  the  negative  portion  of  the  K  range  has  no  effect  on  crack  growth. 
Therefore,  if  the  minimum  K  is  negative,  it  is  assumed  that  AK  equals  the 
maximum  value  of  K.  This  assumption  is  sufficiently  accurate  for  fatigue  life 
prediction  if  the  absolute  value  of  the  calculated  negative  K  is  less  than,  or 
of  the  order  of,  the  maximum  K.  For  autof rettaged  cylinders  operating  at 
relatively  low  pressures  and  for  shallow  cracks,  the  calculated  absolute  value 
of  the  negative  K  at  zero  pressure  is  three  to  five  times  the  maximum  K  at 
pressure.  Such  conditions  would  be  referred  to  as  a  load  ratio,  R  (ratio  of 
minimum  to  maximum  stress  or  K)  of  -3  to  — 5.  Although  there  is  no  known 
fatigue  crack  propagation  data  available  for  such  highly  negative  R  ratios,  an 
examination  of  fatigue  data  on  autof rettaged  cylinders  clearly  indicates  that, 
under  such  conditions,  the  effect  of  the  negative  portion  of  the  K  range  can 
not  be  completely  neglected  in  determining  an  effective  AK. 

As  previously  discussed,  the  maximum  K  is  the  algebraic  sum  of  the  K  due 
to  pressure  (KP)  and  the  K  due  to  the  residual  stresses  (KR) . 

Kmax  =  KP  +  KR 

Since,  for  shallow  cracks  KR  is  negative  and  if  the  negative  portion  of 
the  K  range  is  neglected,  then 
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If  some  fraction  (F)  of  the  negative  portion  of  the  K  range  is  considered 


to  be  part  of  the  K  range,  then 

AK  =  Vx  “  (F)KR  (KR  <  0) 

or 

AK  =  KP  +  KR  -  (F)KR 

=  KP  +  (1  -  F)KR  (7) 


RESULTS 


A  final  equation  for  AK  as  a  function  of  pressure,  residual  stress,  and 
crack  shape  can  now  be  written  by  combining  Eqs.  (4)  through  (7). 


2  24  PY^  Ra  —  RI 

AK  =  (SF)/ia(-~ - +  (SB)  [1 .12  -  -68( — - )  ]  ( 1-F)  >  (8) 

Yz  -  1  RZ  -  RI 


where 

SB  =  -0.4(SY)  +  0.5[otRI  +  0 .4(  SY)  ]  (3) 

The  total  fatigue  life  can  be  determined  by  integrating  Eq .  (1)  between 
the  limits  of  the  initial  crack  depth,  aQ,  and  a  final  value  of  "a"  at  which  K 
equals  a  critical  value.  This  could  be  accomplished  by  substituting  AK  from 
Eq.  (8)  into  Eq .  (1)  and  performing  a  closed  form  integration  between  the 
appropriate  limits.  Since  Ra  in  Eq .  (8)  is  a  function  of  "a",  the  closed  form 
integral  becomes  a  very  complicated  equation.  However,  if  it  is  assumed  that 
Ra  is  constant  over  some  small  interval  of  crack  growth,  and  that  the  constant 
"m"  in  Eq.  (1)  =  3.0,  AN,  the  number  of  cycles  required  to  grow  the  crack  a 
distance  Aa,  or  (a^  -  a(i+l))  can  readily  determined  by 


1  1 

2( - ““> 

a  i+1  da  ^i+1) 

AN  =  /  - - -  (9) 

ai  C(/iFy)3  C( Fy) 3 
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where  Fy  is  defined  as  AK//a  as  determined  from  Eq.  (8).  This  calculation  can 
be  repeated  for  successive  intervals  of  Aa  using  the  average  value  of  Ra  for 
that  specific  interval.  By  summing  the  values  of  AN  for  all  intervals  of  "a" 
from  aQ  to  the  critical  "a”,  the  total  fatigue  life  can  be  determined.  The 
results  of  this  calculation  will  converge  to  the  value  that  would  be  obtained 
by  a  closed  form  integration  as  the  value  of  Aa  approaches  zero.  A 
computational  scheme,  which  can  easily  be  accomplished  on  a  desk  top  or 
personal  computer,  is  outlined  as  follows. 

1.  Using  Eq.  (2)  and  an  iteration  process  to  minimize  determine  the 

value  of  R  at  which  =  0.  This  is  RZ.  Also,  in  the  same  process, 

calculate  otR  at  R  =  RI,  (atRi)  (this  will  be  a  negative  number). 

2.  If  ^  ~0.4,  then  SB  =  °tRI • 

3.  If  otRj_  <  -0.4,  then  calculate  SB  from  Eq .  (3). 

4.  Define  Aa  as  the  cylinder  wall  thickness  divided  by  some  number  of 
intervals . 

5.  Select  an  initial  crack  depth,  aQ,  and  define  the  first  value  of  R  as 
RI  +  aQ  +  Aa/2. 

6.  Calculate  AK  from  Eq.  (8)  using  the  value  of  R  determined  above  as 
Ra.  A  value  of  (SF)  of  0.7  is  suggested  and  the  value  of  F  to  be  used  will  be 
discussed  later. 

7.  Calculate  a  =  Ra  -  RI  and  Fy  =  AK/ /a. 

8.  Calculate  a-^  -  a  -  Aa/2  and  a(x+x)  =  a  +  Aa/2. 

9.  Calculate  AN  from  Eq .  (9).  For  values  of  ”m"  other  than  3.0  a 
different  form  of  Eq .  (9)  is  required. 
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10.  Let  Ra  =  +  ^a  an(^  repeat  steps  6  and  7,  summing  AN  values. 

11.  Continue  incrementing  Ra  by  Aa  and  summing  ANfs  until  AK  exceeds  KC 
for  the  material.  KC  may  be  the  plane-strain  fracture  toughness,  K-j-c,  or  a 
value  somewhat  larger.  Actually,  Kmax  and  not  AK  should  be  used  for  this 
comparison,  but  since  the  crack  growth  rate  is  so  large  at  this  point,  these 
details  have  very  little  effect  on  the  total  life  calculated. 

12.  The  value  of  ^AN  determined  above  is  the  calculated  fatigue  life. 

NUMERICAL  RESULTS 

A  computer  program  to  perform  the  above  calculations,  has  been  written  in 
Basic  for  use  on  an  "Apple”  computer.  A  copy  of  the  program  is  given  in  the 
Appendix . 

In  order  to  determine  the  accuracy  with  which  the  proposed  method  can 
predict  fatigue  lives,  particularly  the  effects  of  diameter  ratio  and 
autof rettage,  a  comparison  is  made  with  the  experimental  results  of  Davidson, 
et  al  (ref  16).  This  is  shown  in  Figure  2.  This  reference  was  selected  for 
comparison  because  it  includes  the  widest  range  of  the  variables  of  interest. 
The  fatigue  lives  for  both  autof rettaged  and  non-autof rettaged  cylinders  were 
calculated  using  the  following  constants  and  material  properties. 

SY  =  160  Ksi  SF  =  0.7 

KC  =  200  Ksi/In  F  *»  0.3 

a0  =  0.005  in.  C  =  3.4  x  10"10 

ra  =  3.0 

■^Davidson,  T.  E.,  Eisenstadt,  R. ,  and  Reiner,  A.  N. ,  "Fatigue  Characteristics 
of  Open-End,  Thick-Walled  Cylinders  Under  Cyclic  Internal  Pressure,”  ASME 
Paper  No.  62-WA-164,  J.  Basic  Engineering,  Series  D,  Vol.  85,  1963,  p.  555. 
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The  crack  depth  interval  over  which  AK  is  assumed  to  be  constant  is 
defined  as  the  wall  thickness  divided  by  some  number,  M.  To  test  the 
convergence  of  the  numerical  integration  technique,  several  cases  were  run 
with  M  varying  from  10  to  100.  A  plot  of  these  results  indicates  that  the 
value  of  life  at  M  =  40  is  less  than  0.5  percent  less  than  the  value  at  M  = 
100.  Therefore,  M  =  40  was  used  for  the  calculations  shown  in  Figure  2. 

The  selection  of  the  values  aQ  =  0.005  and  F  =  0.3  is  somewhat  arbitrary 
and  was  admittedly  made  to  obtain  agreement  with  the  experimental  results. 
However,  the  value  of  aQ  is  a  reasonable  estimate  of  the  maximum  material 
defect  size  that  would  be  expected  in  the  material  used.  The  arbitrary 
selection  of  these  values  would  not  influence  the  agreement  seen  for  the 
slopes  of  the  curves  or  for  the  trend  of  diameter  ratio  effect  on  the 
autof rettaged  cylinders.  The  selection  of  F  does  influence  the  difference 
between  the  autof rettaged  and  non-autof rettaged  cylinders.  However,  for 
reasonable  changes  in  the  value  of  F  from  0.1  to  0.5,  the  resulting  change  in 
stress  at  a  given  calculated  life  is  only  about  15  percent  for  the  2.0 
diameter  ratio,  and  less  for  smaller  diameter  ratios.  This  is  very  much  less 
than  the  difference  between  the  autof ret taged  and  non-autof rettaged  results. 

The  results  of  the  proposed  method  were  also  compared  with  the  results  of 
Throop  and  Fujczak  (ref  17)  which  are  the  only  known  experimental  results 
which  systematically  studied  the  effects  of  crack  shape  and  extent  of 
autof rettage  (overstrain),  and  obtained  fatigue  crack  growth  data.  This  work 

•^Throop,  J.  F.  and  Fujczak,  R.  R. ,  "Strain  Behavior  of  Pressurized  Cracked 
Thick-Walled  Cylinders,"  Experimental  Mechanics,  Vol.  22,  August  1982, 
pp.  277-286. 
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utilized  fairly  deeply  notched  cylinders  which  forced  the  fatigue  cracks  to 
grow  with  specific  crack  shapes.  For  crack  shape  factors  which  would  be 
appropriate  for  these  specific  shapes,  according  to  the  various  references 
previously  discussed,  the  calculated  crack  growth  rates  are  significantly 
greater  than  those  measured.  This  was  noted  by  Parker,  et  al  (ref  13)  and 
others,  and  some  possible  explanations  were  discussed. 

In  order  to  obtain  better  agreement,  shape  factors  for  the  three  crack 
shapes  used  by  Throop  and  Fujczak  were  obtained  by  using  the  proposed  computer 
program  and  adjusting  the  shape  factors  until  the  calculated  fatigue  lives  for 
the  non-autof ret taged  cylinders  agreed  with  the  experimental  fatigue  lives. 
This  procedure  gave  shape  factors  of  0.65  for  the  long  crack,  0.52  for  the 
semi-elliptical  crack,  and  0.37  for  the  semi-circular  crack.  These  values 
were  used  to  calculate  the  crack  growth  curves  for  the  autof rettaged  cylinders 
and  the  results  are  shown  in  Figure  3  along  with  the  experimental  results  of 
Throop  and  Fujczak.  The  proposed  life  prediction  method  generally  under¬ 
estimates  the  life  of  the  notched  cylinders  by  10  to  30  percent.  This  is 
considered  good  agreement  considering  the  inherent  variability  of  this  type  of 
data.  For  the  60  percent  overstrain,  semi-circular  notch  specimen,  the 
prediction  slightly  overestimated  the  life.  However,  there  is  reason  to 
suspect  that  this  experimental  result  may  be  lower  than  expected  due  to  an 
additional  defect  in  the  cylinder. 


■^Parker,  A.  P. ,  Underwood,  J.  H. ,  Throop,  J.  F.,  and  Andrasic,  C.  P. ,  "Stress 
Intensity  and  Fatigue  Crack  Growth  in  a  Pressurized,  Autof rettaged  Thick 
Cylinder,"  ASTM  STP  791,  pp.  I-216-I-237. 
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CONCLUSIONS 


The  fatigue  life  of  a  thick-walled  cylinder  subjected  to  a  constant 
amplitude,  cyclic,  internal  pressure  can  be  calculated  with  reasonable 
accuracy  by  means  of  a  relatively  simple,  fracture  mechanics  based  analysis. 
This  analysis  considers  the  effects  of  crack  shape,  autofrettage  residual 
stresses,  corrected  for  the  Bauschinger  effect,  and  the  negative  values  of  the 
minimum  stress  intensity  factor  produced  by  the  residual  stresses. 

The  computations  required  for  this  analysis  can  be  performed  on  a  micro¬ 
computer,  and  a  "Basic"  program  for  performing  these  computations  is  provided. 
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APPENDIX 


BASIC  ("APPLESOFT")  Program  to  Calculate  Fatigue  Life 
of  a  Pressurized  Cylinder 


10  PRINT  "SELECT  REENTRY  POINT" 

11  PRINT  " 1 =DI AMETER  RATIO" 

12  PRINT  "2=YIELD  STRENGTH" 

13  PRINT  "3=0VERSTRAIN" 

14  PRINT  "4= INITIAL  CRACK  DEPTH" 

15  PRINT  "5= INNER  RADIUS" 

16  PRINT  "6=PRESSURE" 

17  PRINT  "7=END" 

20  INPUT  RE 

23  PRINT  "IF  VALUE  OK  TYPE  0  OR  SPACE, 

50  PRINT  "DIAMETER  RATIO  =  ";W 

51  INPUT  X:  IF  X  THEN  W  =  X 

70  PRINT  "YIELD  STRENGTH  IN  KSI  =  " ; SY 

71  INPUT  X:  IF  X  THEN  SY  =  X 

80  PRINT  "INNER  RADIUS  IN  INCHES  = 

81  INPUT  X:  IF  X  THEN  RI  =  X 
100  INPUT  "OVERSTRAIN  RATIO  =  " ; OS 


IF  NOT  ENTER  NEW  VALUE" 


;Ri 


110 

RZ  = 

i 

1  20 

WE  = 

W  /  (  (W  - 

1  ) 

■* 

□s 

+ 

i ) 

1 30 

WP  = 

W  /  WE 

140 

LWP  = 

=  LOG  (WP) 

150 

R  = 

(WP  +  1)  / 

•n, 

160 

R5  = 

(WE  2  + 

1 ) 

/ 

(2 

* 

WE 

2) 

170 

R6  = 

(1  +  (W  / 

R) 

A 

2) 

/ 

(W 

175 

R7  = 

(WE  2  - 

1) 

/ 

(2 

* 

WE 

A  2) ) 

180 

LR  = 

LOG  (R) 

190 

RS  = 

R5  -  (LWP 

- 

LR) 

- 

R6 

* 

R7 

-  1  ) 

+  LWP 


210  IF  R  =  1  GOTO  300 
220  RX  =  ABS  (RS) 

230  IF  RX  <  .00001  THEN 
235  R  =  R  —  (R  #  . 3  *  RS) 
237  GOTO  160 


GOTO  240 


240  RZ  =  R 

279  PRINT 

290  R  =  l:  GOTO  160 

300  SB  =  RS 

302  IF  SB  >  -  0.4  GOTO  305 

304  SB  =  -  0.4  +  0.5  *  (SB  +  0.4) 

305  PRINT  "RES  BORE  S/SY  R/RI  AT  RES  S=0" 

310  PRINT  SB, RZ 

311  PRINT 

312  SB  =  SB  *  SY 

313  RZ  =  RZ  *  RI 

314  R2  =  RI  *  W 
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320 

321 

335 

340 

341 

o5U 

351 

379 

380 

381 

382 

385 

386 

390 

395 

400 

406 

410 

415 

416 

417 

418 

420 

430 

440 

442 

444 

450 

451 

452 

453 

460 

465 

466 

475 

480 

485 

490 

495 

499 

500 

505 

506 

510 

650 


PRINT  "INITIAL  CRACK  DEPTH  =  " ; AO 

INPUT  X:  IF  X  THEN  AO  =  X 

INPUT  "TO  NOT  PRINT  K  VS  R  TYPE  1 " ; PK 

PRINT  "NO.  OF  INTEGRATION  INTERVALS  =  " ; M 

INPUT  X:  IF  X  THEN  M  =  X 

PRINT  "CRITICAL  K  IN  KSI  IN'- 1/2  =  " ;  KC 

INPUT  X:  IF  X  THEN  KC  =  X 

PRINT  "PRESSURE  IN  KSI  =  ";P 

INPUT  X:  IF  X  THEN  P  =  X 

PRINT  "SHAPE  FACTOR  =  " ; SF 

INPUT  X:  IF  X  THEN  SF  =  X 

PRINT  "RESIDUAL  K  REDUCTION  FACTOR  =  " ; FK 


INPUT  x: 

IF 

X  THEN 

FK  =  X 

SA 

=  (2  * 

P  *  W  A 

2) 

/  (W 

TN 

=  0 

DR 

=  RI  * 

(W 

-  1) 

/ 

M 

R  = 

RI  + 

(DR 

/  2) 

+ 

AO 

IF  RZ  =  RI  THEN  RZ  =  RZ  +  1 
PRINT  CHR$  (9) "BON" 

IF  FK  =  1  GOTO  418 

PRINT  "CRACK  DEPTH  K  PRESSURE  K  RESIDUAL  K  TOTAL 

N  CYCLES" 

XC  =  (2.24  *  F'  *  W  2)  /  (W  "•  2-1) 

XA  =  1.12  -  .68  *  <R  -  RI)  /  (RZ  -  RI) 

XB  =  3. 142  *  (R  -  RI ) 

KP  =  XC  *  (  SQR  (XB))  *  SF 
KR  =  SB  *  XA  *  (  SQR  (XB))  *  SF 
K  =  KP  +  KR  *  (1  -  FK) 

FY  =  SF  *  1.7725  *  (XC  +  SB  *  XA  *  (1  -  FK) ) 

A  =  R  -  RI 

A1  =  A  -  DR  /  2:A2  =  A  +  DR  /  2 

DN  =  ((1  /  SQR  (Al)  -  1  /  SQR  (A2)  )  *  10  10)  /  (1.7  *  FY  ^  3) 

TN  =  TN  +  DN 
IF  PK  =  1  GOTO  475 
PRINT  A, KP, KR, K, TN 
R  =  R  +  DR 
IF  K  >  KC  GOTO  495 
IF  R  >  R2  GOTO  495 
GOTO  420 

PRINT  CHR$  (9) "I" 

PRINT  :  PRINT  "N  CYCLES  BORE  STRESS  INT" 

PRINT  TN, SA:  PRINT 

PRINT  "REENTRY  POINT  =  (7  TO  END) " ; RE 
INPUT  X:  IF  X  THEN  RE  =  X 
ON  RE  GOTO  50 ,70,1 00 , 320 , 80 , 379 , 650 
END 
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TECHNICAL  REPORT  INTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 


CHIEF,  DEVELOPMENT  ENGINEERING  BRANCH 

ATTN:  DRSMC-LCB-D  1 

-DP  1 

-DR  1 

-DS  (SYSTEMS)  1 

-DS  (ICAS  GROUP)  1 

-DC  1 

CHIEF,  ENGINEERING  SUPPORT  BRANCH 

ATTN:  DRSMC-LCB-S  1 

-SE  1 

CHIEF,  RESEARCH  BRANCH 

ATTN:  DRSMC-LCB-R  2 

-R  (ELLEN  FOGARTY)  1 

-RA  L 

-RM  2 

-RP  1 

-RT  1 

TECHNICAL  LIBRARY  5 

ATTN:  DRSMC-LCB-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  UNIT  2 

ATTN:  DRSMC-LCB-TL 

DIRECTOR,  OPERATIONS  DIRECTORATE  1 

DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  L 


NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  WEAPONS  LABORATORY,  ATTN:  DRSMC-LCB-TL, 

OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARMY 
RESEARCH  &  DEVELOPMENT 

ATTN:  DEP  FOR  SCI  &  TECH  1 

THE  PENTAGON 
WASHINGTON,  D.C.  20315 

COMMANDER 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN:  DTIC-DDA  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 

COMMANDER 

US  ARMY  MAT  DEV  &  READ  COMD 
ATTN:  DRCDE-SG  1 

5001  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333 


COMMANDER 

ARMAMENT  RES  &  DEV  CTR 
US  ARMY  AMCCOM 

ATTN:  DRSMC-LC(D)  1 

DRSMC-LCE(D)  1 

DRSMC-LCM(D)  (BLDG  321)  1 

DRSMC-LCS(D)  1 

DRSMC-LCU(D)  1 

DRSMC-LCW(D)  1 

DRSMC-SCM-0  (PLASTICS  TECH  1 

EVAL  CTR, 

BLDG.  351N) 

DRSMC-TSS(D)  (STINFO)  2 

DOVER,  NJ  07801 


DIRECTOR 

BALLISTICS  RESEARCH  LABORATORY 
ARMAMENT  RESEARCH  &  DEV  CTR 
US  ARMY  AMCCOM  1 

ATTN:  DRSMC-TSB-S  (STINFO) 

ABERDEEN  PROVING  GROUND,  MD  21005 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 

ATTN :  DRSXY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005 


NO.  OF 
COPIES 

COMMANDER 
US  ARMY  AMCCOM 

ATTN:  DRSMC-LEP-L(R)  1 

ROCK  ISLAND,  IL  61299 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  (MAT  SCI  DIV)  1 

ROCK  ISLAND,  IL  61299 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACTV 
ATTN:  DRXIB-M  1 

ROCK  ISLAND,  IL  61299 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMD  1 

ATTN:  TECH  LIB  -  DRSTA-TSL 

WARREN,  MI  48090 

COMMANDER 

US  ARMY  TANK-AUTMV  COMD  1 

ATTN:  DRSTA-RC 

WARREN,  MI  48090 

COMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN ,  MECH  ENGR  DEPT  1 

WEST  POINT,  NY  10996 

US  ARMY  MISSILE  COMD 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 

REDSTONE  ARSENAL,  AL  35898 

COMMANDER 

US  ARMY  FGN  SCIENCE  &  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E . 

CHARLOTTESVILLE,  VA  22901 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 

US  ARMY  .AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  DRSMC-LCB-TL , 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT’D) 


NO.  OF 
COPIES 


COMMANDER 

US  ARMY  MATERIALS  &  MECHANICS 
RESEARCH  CENTER 
ATTN:  TECH  LIB  -  DRXMR-PL 

WATERTOWN,  MA  01272 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 
ATTN:  CHIEF,  IPO 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB 

2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CTR 
ATTN:  TECHNICAL  LIBRARY 

CODE  X212 

DAI1LGREN,  VA  22448 


DIRECTOR 

US  NAVAL  RESEARCH  LAB 
2  ATTN:  DIR,  MECH  DIV 

CODE  26-27,  (DOC  LIB) 
WASHINGTON,  D.C.  20375 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
1  ATTN:  AFATL/DLJ 

afatl/dljg 

EGLIN  AFB,  FL  32542 

METALS  &  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  LAB 
1  505  KING  AVENUE 

COLUMBUS,  OH  43201 


1 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 

US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  DRSMC-LCB-TL , 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


NO.  OF 
COPIES 


1 

1 


1 

1 


1 


